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Abstract 


We show that the radiative corrections to charged current (CC) nuclear reac¬ 
tions with an electron (positron) in the final state are described by a universal 
function. The consistency of our treatment of the radiative corrections with 
the procedure used to extract the value of the axial coupling constant gA is 
discussed. To illustrate we apply our results to (anti)neutrino deuterium dis¬ 
integration and to pp fusion in the sun. The limit of vanishing electron mass 
is considered, and a simple formula valid for Eohs ^ 1 MeV is obtained. 
The size of the nuclear structure-dependent effects is also discussed. Finally, 
we consider CC transitions with an electron (positron) in the initial state and 
discuss some applications to electron capture reactions. 

PACS number(s): 13.10.-|-q, 13.15.-|-g, 25.30.Pt 


I. INTRODUCTION 

The physics of neutrino flavor oscillations has entered a new era. There is now a consensus 
that the observation of the atmospheric neutrinos [1] can be interpreted as an evidence 
that muon neutrinos born in the atmosphere oscillate into tan neutrinos. At the same 
time, numerous measurements of solar electron neutrino fluxes have been pointing toward 
neutrino oscillations ever since the Homestake experiment [2]. The evidence for oscillation 
of solar neutrinos has been strengthened recently with the results from Sudbury Neutrino 
Observatory (SNO) [3,4] that are independent of the solar model (SSM). With the evidence 
for oscillations at hand, the accurate determination of the neutrino masses and flavor mixing 
angles is the most important task for neutrino physics at present. Carrying out the above 
program requires accurate knowledge of the cross sections of the reactions used for neutrino 
detection, and the Standard Model (SM) radiative corrections, which typically shift the 
leading (tree) order values by 3-4%, must be taken into account. 

Complete one-loop SM radiative corrections to the CC and neutral current (NC) deuteron 
disintegration by electron neutrinos used in SNO measurements have been calculated in Refs. 
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[5,6]. We have shown in Ref. [5] that while the differential CC cross section depends on the 
actual detector properties {e.g. on the bremsstrahlung detection threshold the total 

cross section is a detector-independent quantity as long as the hnal state electron is always 
detected, i.e., the total number of the neutrino interactions is determined. In the simple 
case when 0 one hnds: 

dacciEobs) = da'S^S'iEobs) (l + %{Eobs)) , (1) 

where da^^‘^{Eohs) and dacc{Eobs) are, respectively, the leading and the next to leading 
order in a differential cross sections, Eobs is the energy observed in the detector (charged 
lepton energy plus possible bremsstrahlung photon energy), and g{Eobs) is given by the Eq. 
(3) below. It is understood that the cross section da^^^^Eobs) includes the Fermi function 
E{Z, A, Eobs) to account for the distortion of the final state electron wavefunction by the 
Coulomb held of the hnal nucleus (see e.g. Ref. [7]). While a similar approach to the 
treatment of radiative corrections has long been used in the analysis of nuclear beta decay 
[8,9], here we extend it to a wide class of reactions involving neutrinos, and formulate the 
conditions for its applicability Specihcally, we generalize the results of Ref. [5] and analyze 
the following four features of the radiative corrections. 

1. Universality. The function g{E) is universal for a class of nuclear reactions involving 
neutrinos. In particular, it can be used to evaluate the radiative corrections to the total 
cross-sections of reactions that have an electron (positron) in the hnal state. It is applicable 
not only to neutron and nuclear /^-decays, but also to the antineutrino capture reactions, 
low energy CC neutrino-nucleus disintegration, nuclear fusion reactions accompanied by the 
electron(positron) emission, etc. For illustration we consider the following CC processes 
involving two nucleon systems: 


Ue + d^p + p + e , 
z/g + d—, 

p + p^d + e^ + Ue. ( 2 ) 

For reactions where the electron spectrum is narrow (precise conditions are given in Section 
II) we provide a simple and accurate way to calculate the correction to the total cross section 
without the complicated integration over the outgoing electron or positron spectrum. 

In order to demonstrate one of the possible applications of our results, note that the last 
reaction in Eq. (2) is the main reaction for the pp chain that powers the Sun (see e.g. Ref. 
[10]). The total rate of this reaction is, therefore, constrained by the known solar luminosity. 
That rate, together with the cross-section for the pp reaction, is used to calibrate the SSM 
and to predict many other quantities, in particular the flux of the neutrinos studied 
at SNO. Below we show that a proper treatment of the radiative corrections to Spp, the 
S-factor for the pp reaction, lowers the SSM prediction to <I>(®B) by about 0.6% relative to 


^The connection between g{E) and the function G{Em, E) introduced by Sirlin [8] for the descrip¬ 
tion of allowed /3-decays is discussed in Section VI. 
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the currently accepted value [11]. ^ 

2. Axial coupling constant qa- Since the reactions in Eq.(2) are dominated by the nuclear 
axial current, it is important to use a procedure consistent with the dehnition of the axial 
coupling constant qa- We briefly discuss the conventional dehnition of qa based on the 
value extracted from neutron decay and give a prescription for evaluation of the radiative 
corrections consistent with that dehnition. This issue has been discussed in the literature 
(see e.g. Ref. [13]) but we feel that it is sufficiently important to reiterate here. 

3. Nuclear structure-dependent ehects. In deriving universality properties of the radia¬ 
tive corrections we work in the “one-body” approximation (in analogy to the early appli¬ 
cations to the nuclear beta decay). Namely, we evaluate the radiative correction to the 
CC reaction on a single nucleon, and then use it to compute the correction to the reaction 
of interest. Ehectively, only the corrections of the type shown in Fig. 1(a) are included in 
such a calculation. Nuclear structure-dependent many-body corrections of the type shown 
in Fig. 1(b) are thus neglected. We show that a comparison of such corrections to analogous 
ehects contributing to superallowed Fermi transitions in nuclei suggests that they enter at 
0.1% level. As we discuss in Section V, however, this estimate of universality-breaking con¬ 
tributions is not airtight, and a detailed computation is needed to obtain precise numbers. 

4. Collinear singularities. We study the behavior of the radiative corrections in the 
limit rrie —^ 0, where rrie is the electron mass. Separate contributions to the corrections are 
singular in this limit, and until now it has not been explicitly shown in the literature how 
the singularities cancel when all contributions are added. Such cancellation must take place 
according to the well-known Kinoshita-Lee-Nauenberg (KLN) theorem [14], and in Section 
VI we show how this happens. As a corollary we obtain a simple expression for g{E) valid 
for E 3> rrie- In practice, radiative corrections evaluated according to this simple expression 
are accurate to about 0.1% for energies above 1 MeV. 

Although these four topics constitute the main focus of our study, we also discuss pro¬ 
cesses, such as capture reactions, involving charged leptons in the initial state. In this case, 
the radiative corrections do not display the same kind of universality applicable to the re¬ 
actions in Eq. (2). Nevertheless, the treatment of the two cases is sufficiently similar that 
we feel a brief discussion - given at the end of the paper - is warranted. 


II. UNIVERSALITY OF THE RADIATIVE CORRECTIONS TO REACTIONS 
WITH THE ELECTRON/POSITRON IN THE FINAL STATE 

We begin with the hrst reaction in Eq.(2). As was shown in Ref. [5], the total cross 
section is independent of the photon detection threshold In the limit 0 the 

function g{E) in Eq.(l) has two parts: 

g{E) = gy{E) + gb{E) . (3) 


^The “outer” radiative corrections to Spp have been originally evaluated in Ref. [12], whose results 
are, however, not used in Ref. [11]. Results obtained using our simple prescription agree with the 
calculation presented in Ref. [12]. 
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Here, gv{E) (for the virtual part) and gb{E) (for the bremsstrahlung part) are given by the 
following lengthy expressions (the cutoff parameter A that appears in Ref. [5,6] is set equal 
to the proton mass): 


and 
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In the above equations, me is the electron mass and: 


P{E) = 


E^ 


mt 


\E\ 


Jo X 


(4) 


(5) 


( 6 ) 


The closed form for the integrals appearing in Eq. (5) is given in the appendix. Although 
g{Eobs) in general depends on the value of the cutoff parameter A, the choice of A affects 
only the constant, energy independent part of gv{Eobs)- We chose A = Mp and adjusted the 
constant (-0.57) in the Eq.(4) so that the result matches with the corresponding expression 
derived by Sirlin in Ref. [9] based on current algebra (the same approach was chosen in 
Ref. [6]). The dominant uncertainty in g{Eobs) is associated with the value of the matching 
constant and is briefly discussed after Eq. (8). 
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It is remarkable that g{Eobs) is a function of a single parameter - the observed energy 
Eohs - and does not depend on the electron and photon energies separately. Moreover, 
formulas Eq.(4) and Eq.(5) make no reference to the parameters that describe the leading- 
order differential cross section, such as pp scattering length or the effective range (see e.g. 
Ref. [15] for dehnitions). It is, therefore, reasonable to ask whether g{Eobs) is a universal 
function that describes the radiative corrections to a whole class of processes in the limit 
£jmm Below we show that this is indeed the case. 

Consider Erst, in the one-nucleon approximation, contributions from virtual photon ex¬ 
changes. As pointed out in Refs. [5,6], it is convenient to split these contribution into two 
pieces: with k < A and with k > A, where k is the scale for the virtual photon momenta 
and A is a cutoff, taken to be of the order of 1 GeV. Contribntions with k > A, usnally 
called “inner” radiative corrections, are combined with the boson exchange box graphs 
to prodnce an energy-independent contribution to g{Eobs) [6]: 

g^{E,,^Y^^ = 2\n(^^^ . (7) 

This contribntion is obviously common to all reactions in Eq. (2) since it contains no param¬ 
eters that distinguish between them. 

The piece with k < A is more complicated. The virtual exchanges of a low energy photon 
for all three reactions Eq.(2) are shown in Fig.2. The correction to Fig.2(a) is [5,6]: 


gYEobsT^'" = ^\n 


+ 21n 



1 ( l+(3{Eobs) \ ^ 

_2P{Eobs) \^ — P{Eobs)) 


0.57 + SC{A) , 


( 8 ) 


where A is the “photon mass” used as the infrared regulator and A{P) is given in Eq. (4). 
The parameter SC represents a matching constant arising from presently incalculable non- 
perturbative QCD effects. Its A-dependence must cancel the corresponding dependence in 
g{Eobs)- Unfortnnately, no existing calculation of SC has demonstrated this cancellation 
explicitly. In addition, SC will contain A-independent terms that depend on the light-qnark 
masses, Aqcd, etc.. A complete, first principles compntation of these contribntions has yet 
to be performed, so we must rely on models. Thus, in practical calculations we set this 
constant to zero for A = to match the model resnlts of Refs. [9,16,17]. The uncertainty 
in gv{E) associated with SC is estimated in Ref. [17] to be abont 0.08%. 

Even thongh the virtnal momenta in Figs.2(b,c) are different from those on Fig.2(a) their 
contribntions are also given by the Eq.(8). The easiest way to see it is as follows. Consider 
Fig.2(b) first. This graph can be obtained from Fig. 2(a) through time reversal followed by 
replacements e“ (initial) —> e"*" (final) and z/g(final) —z/g(initial). Since time reversal is an 
exact symmetry of QED, the scattering amplitnde is nnchanged np to an overall phase. The 
result of replacements of the particles with antiparticles can be fonnd by ntilizing crossing 
symmetry. However, the replacement pe —^ —pe does not change j3{Ee) (see Eq.(6) for 
the definition). Since Eq. (8) depends on the electron momenta only throngh j3{Ee) it is 
obviously invariant. Hence, Eq. (8) is applicable to the second process in Eq.(2). The 
validity of the Eq. (8) for the third process in Eq.(2) trivially follows from the fact that Eq. 
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(8) is independent of the neutrino 4-nionientuni. Indeed, going from Fig.2(b) to Fig.2(c), 
which is accomplished by replacing the antineutrino with the neutrino and changing the sign 
of its 4-momentum, leaves the Eq. (8) invariant. 

The bremsstrahlung contributions are also identical for all three reactions under consid¬ 
eration [Eq.(2)]. In particular, for these reactions, only the bremsstrahlung from the charged 
lepton is important since all other charged particles are heavy. Therefore, for each of the 
processes under consideration there is only one bremsstrahlung diagram of the type shown 
in Fig.3. The amplitude for this graph has the form: 

M, ~ . ( 9 ) 

where e is the electron charge, is the photon polarization, and and H^{q) represent, 
respectively, the leptonic and the hadronic contributions to the process: 

m=l 

= I ■ ( 10 ) 

Here and are, respectively, the initial and the final states of the nuclear and the 
leptonic parts of the process, is the relevant charged current, is the electromagnetic 
current, and k is the momentum of the emitted photon. The summation index m in the 
hadronic matrix element runs over all nucleons in the nucleus. The dependence of on 
the momentum transfer q is explicitly shown. In the long wavelength approximation, valid 
for low neutrino energies, one can neglect the dependence of on the spatial components 
of q. On the other hand, the dependence of on the time component of q could be 
crucial. Indeed, energy conservation implies that the initial Ei and hnal Ef energies of 
the hadronic system obey the relation Ef = Ei + Clearly, the nuclear matrix element 
depends sensitively on the value of Ef, the excitation energy of the final state. In the case 
of CC neutrino deuteron disintegration, for example, changing Ef causes a dramatic change 
in the shape of the hnal state pp wavefunction and the corresponding change in the size of 
the nuclear matrix element [5]. 

With these consideration in mind, and using the results of Ref. [5,6], it is easy to show 
that the square of the matrix element Eq.(9) averaged over the initial and summed over the 
hnal spins has the form 



{\M^{E^,Ee,E^,x)\‘^) = -{\Mo{E^,Eobs)\‘^) x -^E{E^,Ee,x) , 

^ ^ohs 

q^ = EEu — Eg — E^ = EEy — Eobs ■ (H) 

Here {\M.o{Ey, Eohs)\^) is the spin-averaged square of the matrix element in the leading order 
of perturbation theory and x is the cosine of the angle between the momenta of the photon 
and the electron. The dependence of M.q on q^ is implicit through the second line of Eq.(ll), 
where the sign of the neutrino energy depends on whether the neutrino is in the initial (-I-) 
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or final (-) state. F{E^, Ee,x) is given by the second line of Equation (11) in Ref. [6] 


E{E^,Ee,x) 


2El{E, 


E^ 

- P{E,)Kx) 


(E, + E,\P\E){l-K^xyE^^) 

V 2Ee J {E^ - (3{Ee)Kxf 


( 12 ) 


As before, A is the infrared regulator (“photon mass”). While the form of the matrix element 
Afo may depend on the particular process under consideration, the function E{E^, Ee,x) is 
universal. 

In order to complete the calculation of the bremsstrahlung contribution to the cross 
section it is sufficient to multiply the hrst line of Eq.(ll) by the appropriate phase space 
factor and integrate over the final state momenta to obtain 


da'ljc = ^obs) X 


« fd{Ee)El 


E{E^, Ee, x)dxK 


dKdEe 

dEf^i)^ 


( 13 ) 


where again Eohs = Eg + E^ and da^^‘^{Ey, Eobs) is the leading order differential cross section. 
If there is at least one more particle in the hnal state except the electron and the photon, 
its phase space factor can be used to eliminate the spatial part of the overall momentum- 
conserving delta function, which is always present in the expression for the differential cross 
section. Keeping Eobs fixed, one can then perform integrations over x and Ee. Changing the 
integration variables from {Ee, A'} to {Ee, Eohs} and integrating over Ee from nie to Eobs we 
obtain: 


m. cy P ^obs r 1 

dalciEu,Eobs) = da^S^iE,,Eobs) x - / dEe / dx 

TT J rrie J — 1 

X ^-1 „2 E(Eobs ~ Ee, Ee, x)(Eobs ~ Ee) 

PK^obsji^obs 

= da'p^ciEu: Eobs) X -gbremiEobs) ■ (14) 

TT 

The above integration is delicate due to the presence of infrared divergence. The result, 
however, is well known (see e.g. Ref. [6]): the infrared divergence in gbrem cancels that in 
Eq.(8). 

Now, the final result for the sum of all corrections is 

g{Eobs) = gv{Eobs)'^~^ + gv{Eobs)'^~^ + gbremiEobs) , (15) 

with g{E) given in Eq.(3) for A = Mp. This completes the proof of the universality of this 
function. The dependence of g{E) on energy is shown in Fig.4. 


^Our definition differs from that in Ref. [6] by a factor of 2, which if absorbed in F{E^, Ee, x) to 
simplify the expression in Eq. (11). 
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III. EXAMPLES 


With help of Eq.(l), the function g{Eobs) obtained in the previous section can be used 
to account for the radiative corrections to the differential cross section to all reactions listed 
in Eq.(2). A dramatic simplihcation occurs if one is only interested in computing the total 
cross section. The latter is given by the integral 

^max , 

acc= da^S^{E)(l + %{E)) (16) 

J me \ TT / 


where the subscript “obs” on E is not shown. We argue below that if the leading order 
electron spectrum is sufficiently narrow there is no need to compute the integral explicitly. 

As it is apparent from Eig. 4, g{E) is a smooth fnnction of its argument. If the shape 
of the leading order differential cross section is such that only a certain range of energies 
dominates the above integral, g{E) can be expanded around some point Eq inside that range. 
Snch an expansion leads to the following series for the total cross section; 


ffoc = (l + 

+ -{g\Eo) / {E - E,)da^S^{E) 

TT y Jirie 

+ ^9"{Eo) {E - Eofda^^\E) + ■ ■ ■ j 


‘ the point Eq is chosen in such a way that 


z?. 


„ Tpmax m 

,ff Ed,jJ;fFiE) 


= (E) 


( 17 ) 


(18) 


then the second term in Eq.(17) vanishes. Here, {E) represents the average observed energy 
calculated with the leading order electron spectrnm. The series now has the form 


<^cc = (l + A((-B)) + ^9"((B»(«5t + ■ ■ ■) 
SZ”(E - {E)fdal-S‘(E) 


(19) 


The above formnla shows that the total cross section can be represented as a series 
in moments {SE"') of the leading order electron spectrnm with coefficients given by 
a/{7rn\)g^'^\{E)), with snperscript indicating the nth derivative. Since we expect the fi¬ 
nal uncertainty in the cross section to be of the order of 0.1% (see Section V), the series can 
be truncated at the leading term if the following condition holds: 


rv 

—g"({E)){6E'-‘)<(i.l% 


( 20 ) 


Below we consider some examples for which snch trnncation is possible. 



A. Neutrino and Antinentrino Denterium Disintegrations 


Consider the first two reactions in Eq.(2). The shape of the leading order electron 
spectrum is dominated by the overlap integral of the wavefunctions of the deuteron in the 
initial state and the two nucleon system in the hnal state. The overlap integral depends 
on the relative momentum of the two nucleons, and it falls rapidly when the momentum 
becomes larger than the inverse of the scattering length for the hnal state system (see Ref. 
[5] for discussion). Because of this feature the electron spectrum is strongly peaked near the 
endpoint (see e.g. Ref. [15]), and its width is determined by the corresponding scattering 
length; 




1 


a 


2 

N1N2 


Mp 


( 21 ) 


where aNiN 2 is the scattering length for the hnal state containing nucleons Ni and N 2 . The 
direct evaluation gives SEpp ^ 0.7 MeV, 6Enn ~ 0.2 MeV. 

For Ey > A MeV (as it is the case for the CC reaction at SNO) one can use Eq.(3) to 
show that 


a 

n 


9"{E^) 


0.0037 


( 22 ) 


for Ey is in units of MeV. Therefore, for E^>5 MeV we have for both channels: 

rv 

—g'\{E)){6E^) <0.01% (23) 

Zn 

which is an order of magnitude better then necessary. 

To verify the validity of the truncation procedure we performed a series of calculations 
using Eq.(16) for the neutrino energies E^ > 4 MeV for both pp and nn hnal states. The 
truncated expression for the total cross-section 

<€ = (1 + (24) 

was valid to within the estimated error Eq.(23). For instance, for the pp hnal state at 
E^ = 4 MeV the exact one-loop radiative correction to the total cross section is equal to 
3.28%, whereas Eq.(24) gives 3.27% with {E) = 2.35 MeV. 


B. The pp fusion reaction 


The electron spectrum for this reaction [the third in Eq. (2)] is essentially determined 
by the hnal state phase space factor: 

-^^P{EM{Eo-E^f (25) 

where Eq is the maximum electron energy. Such shape is also typical in nuclear beta decays. 
For the pp reaction in the Sun, the kinetic energy of the initial state protons may be neglected 
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in computing Eq. In that case, Eq ^ 2Mp — Md = 0.93 MeV. For the spectrum given by the 
above equation, the first two moments are well approximated by 


(E) 

{6E^) 


Eq + rrie 


{Eq - rUef 
24 


(26) 


Also, one can show from Eq.(3) that for E^ < 1 MeV one has to a good approximation: 


a 0.0031 

-g'\E,) ^ - -^ 

71 me{Ee - rrie) 


(27) 


with Ep and rrip in units of MeV. Therefore 


a 


g"{{E)){5E‘^) ^ 0.06 


{Eq - me), 


mp 


(28) 


In particular, for pp reaction we get the truncation error estimate of 0.05%, which is more 
than satisfactory. 

The rate of the pp reaction in the Sun strongly affects the flux of neutrinos ^{^B) 
recently detected at SNO [3,4] The relationship between ^(^B) and the S-factor for the pp 
reaction is [11]: 


4>(®B) 


rsj 


C-2.6 

^PP 


( 29 ) 


Since the fractional change in the total pp cross section due to the radiative corrections 
translates into the same change in Spp we obtain: 


5 $ ( 85 ) 6 Spp 6aPP 

—f —L = _2.6—^ = -2.6- 

S„ a-- 


(30) 


The exact evaluation using Eq. (16) gives 3.87% while the first term in Eq.(19) gives 3.86% 
for the average electron energy of 0.67 MeV. We see that the truncation error is well within 
the aforementioned 0.05% estimate. 

To compare our calculation with the result adopted in Ref. [11] we need to subtract the 
“inner” part of the correction to isolate the “outer” piece: 


g({E)r“'= s({E)) - f 


( 31 ) 


To emphasize that the inner correction is independent of the electron energy we do not write 
the argument for We follow the convention adopted in Ref. [13] and identify the outer 

piece with the contribution from Sirlin’s function G{Em, E) (see Ref. [8]) integrated over the 
electron spectrum. To isolate the remaining inner piece it is convenient to use Eq. (29) of 
Ref. [6]. Dropping G{Em,E) on the RHS of Eq. (29) we obtain: 

pinner _ 2 In + 0.55 (32) 
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where the 0.55 on the RHS is obtained in Ref. [6]. R includes non-asymptotic contributions 
from the weak axial vector current (denoted by C in Ref. [6]) as well as perturbative QCD 
corrections. The value of the inner radiative correction is evaluated to be 2.25%. Since it is 
independent of the electron energy we can simply subtract it from the total correction to hud 
the outer piece. The latter is, therefore, equal to 1.62%, which is 0.22% higher than 1.4% used 
in Ref. [11]. The difference is likely caused by a slightly lower value of Eq in the pp reaction 
compared to the corresponding quantity in neutron decay: = Mn — Mp = 1.3 MeV 

compared to 0.93 MeV for the pp. Since g{E) is a decreasing function of E, smaller Eq in 
the pp reaction means a larger correction. Using the actual value of the radiative correction 
reduces the prediction for by 0.57% relative to the present value. Although not very 

signihcant, this shift is comparable to the total theoretical uncertainty of 0.5% assigned to 
Spp in Ref. [11]. 


IV. AXIAL CURRENT COUPLING CONSTANT 

An accurate prediction of the cross section or rate of a process that depends on the nuclear 
matrix elements of the weak axial vector current requires an accurate knowledge of the axial 
coupling constant qa of the nucleon. If radiative corrections are included one must use 
a renormalization scheme consistent with the procedure used to extract the experimental 
value of this constant. Here we briefly review the prescription for evaluating radiative 
corrections to the hadronic matrix elements of the axial current when using the value for qa 
recommended by the Particle Data Group [18]. 

At leading order in electroweak perturbation theory (but to all orders in the strong 
interaction) the nucleon matrix element of the charge-raising current has the form: 

{Nf,Pf\^qklt,{^-'lb)T^qk\Ni,Pi) = Pigvl^ + 

O 

-S'aTmTs - gpq^ilb)n , (33) 

where the gis represent the corresponding coupling constants, and the sum on the left hand 
side runs over all quark flavors. The relevant form-factors are evaluated at q^ = {pf—PiY = 0) 
and only contributions from the first class currents are included. The above formula leads 
to the inverse neutron lifetime 

% ^ ~ ^ud^’p (dv + H-) (34) 

where Vud is the first element of the CKM quark flavor mixing matrix, is the Fermi 
coupling constant determined through the muon lifetime, and ellipses represent the small 
pseudoscalar and weak magnetism corrections. In the absence of isospin breaking the con¬ 
served vector current hypothesis (CVC) holds, and gy = 1. We will not discuss corrections 
to this approach in the following (see e.g. Ref. [19]). On the other hand, the value of gA is 
not protected from renormalization by QCD effects. For the purpose of this discussion we 
will call gA appearing in Eq.(33) the “fundamental” axial coupling and accentuate it by a 
circle. 
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In practice, one uses the combination of measurements of the neutron lifetime and its 

O 

parity-violating decay asymmetry, which depends on the ratio fl'A/fl'y- If CVC is invoked, one 
can also determine Vud- However, for the accurate determination of all involved parameters 
the above formula for the neutron lifetime is insufficient and the electroweak radiative cor¬ 
rections must be taken into account. The radiative corrections to neutron decay have been 
evaluated in a number of papers (see, e.g., Ref. [9,13,20]). The general result of these papers 
is that the neutron differential decay rate, averaged over the neutron spin and integrated 
over all variables except the electron energy, has the form 


dTn 

dEe 


GiV 


F''ud 


[1 -|- 2Av{Ee)] gY+[l + 2A^(i?e)] ^9^^ 


(35) 


where Eg is the energy of the emitted electron, and Av{Ee) and Aa^E^) are, respectively, the 
radiative corrections to the vector and the axial vector current contributions to the decay. It 
is important to note that these corrections are in general not equal. However, for the level of 
precision we are interested in their difference can be regarded as a constant independent of 
the electron energy [20]. Although the value of this constant can in principle be computed 
in Standard Model, the practical calculation contains hadronic structure uncertainties [5]. 
It is, however, possible to absorb such uncertainties in the modihed dehnition of gA- 


dVn 

dEe 


~ [1 -h 2Av{Ee)] G%V^^ {gy + 


1 


gA = fi'A [1 + AA^Ee) — Av{Ee)] — 9 a{^ + d) . 


(36) 


The advantage of the above dehnition is that the same radiative correction Ay{Ee) domi¬ 
nates the corrections to the rate of superallowed nuclear beta decays, which are pure Fermi 
(vector current) transitions. Its uncertainty contributes to the uncertainty in the value of 
Vud extracted from such decays. At present, the most conservative estimate of the combined 
experimental and theoretical uncertainty in Vud is 0.07% [18], which sets the upper bound on 
the uncertainty in Av{Ee). With Av{Ee) constrained in this way, gA dehned in the second 
line of the above equation can be experimentally determined from the neutron lifetime and 
its decay asymmetry. It is precisely the value for gA dehned in the above equation, not for 

O 

the fundamental 9 a, that is quoted in the Particle Data Group [18]. The current number is 
gA = 1.2670 ± 0.0030. 

The above considerations -in particular Eq.(36) - show that instead of evaluating the 
radiative corrections to the hadronic matrix elements of the axial current one can use the 
corrections to the corresponding matrix elements of the vector current in combination with 
the modihed axial coupling constant gA- For the reactions in Eq.(2), there is no need to 

O . 

know the fundamental coupling Such an approach eliminates most of the hadronic un¬ 
certainties mentioned in Ref. [5], where the discussion was given in terms of the fundamental 

O 

9a. 


^For other processes, however, one does require 9a- In neutral current lepton-nucleon scattering, 
for example, one must include only radiative corrections to the neutral current amplitude. The use 

O 

of gA rather than Qa would erroneously introduce the effects of charged current radiative corrections 
to the neutral current amplitude. 
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V. NUCLEAR STRUCTURE-DEPENDENT CONTRIBUTIONS 


Up to this point, our discussion of radiative corrections has applied to the one-nucleon 
approximation; corrections were computed assuming that both the W boson and the photon 
couple to the same nucleon. This approximation corresponds to graphs in Fig. 1(a). There 
are, however, other contributions, such as those in Fig. 1(b). These contributions depend 
on the nuclear structure, which, in general, renders them non-universal. Below we present 
arguments suggesting that such corrections should contribute to the cross section at the 
level of 0.1%. 

If the nuclear structure-dependent contributions to the cross section are not neglected 
then Eq.(l) is modihed to [16]: 

dacciEobs) = da^S'iEobs) (l + ^ [giEobs) + C^s]) , (37) 

where the quantity Cns represents corrections that depend on nuclear structure. Its ac¬ 
curate evaluation requires knowledge of the full nuclear propagator between the weak and 
electromagnetic current insertions. Although possible in principle, the calculation of this 
quantity is a complex problem even for such a simple nucleus as the deuteron. Before dis¬ 
cussing Cns for the deuteron it is instructive to turn to the case of superallowed nuclear 
Fermi (3 decays where the analogous term has already been studied. 

Nuclear structure-dependent radiative corrections to superallowed nuclear f3 decays can 
be split into two categories: those induced by weak vector current (VC) and by weak axial 
vector current (AC). For the pure Fermi transitions, the contributions from the VC are 
suppressed by powers of the electron energy or first generation quark mass [21] and are 
negligible. In Refs. [16,17] contributions to Cjjg induced by the AC have been evaluated for 
Fermi transitions in a number of nuclei using a shell model approach. was modeled by 
contributions analogous to those shown in Fig. 1(b). According to Table II of Ref. [17] the 
magnitude of Cjjg never exceeds 1.348. 

The results presented in Ref. [16,17] are, of course, model-dependent. However, for pure 
Fermi transitions within the Standard Model there exist indications that such approach 
gives fairly reliable values for Cjjg. Indeed, these values (along with some other nucleus- 
dependent corrections) successfully bring /t-values of various superallowed nuclear beta 
decays in agreement with each other (at about the 0.1% level) as required by CVC (see 
e.g. Ref. [17] for discussion). This success represents a non-trivial test of theoretical nuclear 
structure corrections in superallowed nuclear beta decays and suggests that even if the 
model calculations of such corrections differ from their actual values the discrepancy must 
be a nucleus-independent systematic effect. Although existence of such an effect might help 
to solve the CKM unitarity problem within the Standard Model, the authors of Ref. [17] 
argue that it does not appear plausible under any reasonable circumstances. 

These considerations suggest that the existing model calculations for superallowed decays 
produce reliable values for for a large number of nuclei with 10 < A < 74 [17]. Therefore, 
one might expect that application of similar techniques to the transitions in Eq. (2) will also 
produce reliable results. Up to terms suppressed by electron energy or light quark masses, 
the dominant nuclear structure-dependent effects in this case arise from the VC, rather 
than the AC as for Fermi transitions [21]. Despite this difference, one might nevertheless 
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argue that an analogous model calculation of for the deuteron should be comparable 
in magnitude to up to appropriate scale factors (see below). According to Ref. [16] 

magnitude of is related to the magnitude of the typical velocity of a bound nucleon: 


(pjv) _ M 

^NS ~ “ 


M, 


N 


C 


(38) 


where Mjv is the nucleon mass, is the characteristic nucleon momentum {e.g. Fermi 
momentum), and is the corresponding nucleon velocity . This utv- dependence follows 
directly from the expressions for the nucleon weak and electromagnetic currents in the non- 
relativistic case. 

A similar situation occurs for Gamow-Teller transitions. The dominant nuclear structure- 
dependent contribution arises from the antisymmetric correlator of hadronic vector currents 
appearing in one-loop amplitude [21]: 

VV = i /<iVe-“"’(/|r{U(i)y™(0)}|i) , (39) 


where he and Uy are Dirac spinors for the electron and the neutrino, Vx and J™ are, 
respectively the nuclear weak vector current and the electromagnetic current operators, and 
/ and i represent the hnal and the initial nuclear states. Since two indices of the 
tensor are contracted with the indices of the two hadronic vector currents at least one of 
these indices must be space-like. Moreover, in the non-relativistic limit relevant for low- 
energy nuclear dynamics, the spatial components of nucleon vector currents are 0{{vn)/ c). 
Hence, we would also expect to scale with v^. 

Since the characteristic nucleon momentum in the deuteron is signihcantly lower than 
that in all nuclei listed in Table II of Ref. [17] our naive scaling argument would give a 
smaller magnitude for G^|((i): G^|((i)/G^|(A) Pd/PA- Here, the quantities describing 
the deuteron and the nucleus are accompanied by d and A, respectively. For atomic number 
A ^ 1 the Fermi momentum is about 300 MeV [7]. On the other hand, in the deuteron the 
corresponding momentum scale is simply given by its inverse size, 7=45 MeV [15]. Taking 
the largest value for from Table H in Ref. [17] and rescaling it by the corresponding 
momentum ratio we get the following estimate for the maximal expected size of for the 
deuteron disintegration: 

~ ^ X ^ X 1-348 ^ 0.05% (40) 

Although such a correction is negligible at present, we cannot rule out fortuitous few-body 
structure effects that might lead to an enhancement. Completion of a detailed few-body cal¬ 
culation would determine whether any such enhancement occurs and provide hrmer bounds 
on the theoretical, nuclear structure-dependent uncertainty. 


VI. TAKING THE ELECTRON MASS TO ZERO 

According to KLN theorm [14], a properly dehned scattering cross section should have 
no external mass singularity. In this section we demonstrate how the result Eq.(3) is in 
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agreement with this requirement. We take rrie —^ 0 and show that the function g{E) in 
Eq.(3) is hnite in this limit. The results of this section also have a practical application 
since we obtain a simple expression for g{E) valid for E 3> rrie, which allows one to evaluate 
radiative corrections with accuracy of about 0.1% for energies above 1 MeV. 

In the limit me —0 we have: 


13(E) 

1 + /3 

m 

L{-P) 


1 - ml/2E^ 

. 2E 

2 In — 

me 

T A + /^\ 


TT 

12 

0 


' 2/3 ' 

1 + /3 


TT 

’Y 


(41) 


Using these expression it is straightforward to take the limit me —0 in Eqs. (4) and (5). 
We have 


E 


g^{E) ^ 21n 


me 

Me, 


E 


TT 


A{I3) ^ - In^-2 In 2 - - In — + In 2(1 - In 2)-1 


+ xln 


TT 


2 

fMp 
\me 


2 1 


m, 

+ In 


6 




me 


me 


ln2(l - ln2)-0.57 

6 4 


(42) 


The integrals in Eq.(5) can now be evaluated analytically: 


2E^p{E) 


[^{E-x) In ( 

Jme \ 


1 + P{x)' 
1- p{x)^ 


dx 


^ (In—+ ln2-? 


E(3{E) 


xp{x)E{x) - Ep{E)E{E) 
rUe E - X 


dx 


me 

E TT^ 

-2 In-21n2-+ 4 

me 3 


(43) 


Together with the limit for C{j3) 


E 


E 


TT 


C(/3) ^ -In^ — + ln-ln2(l -ln2)-+ 1 

me me 6 


(44) 


we have 


, , <2 E \ E 

gb{E) ^ - In-- In — 

me 2 me 

-ln2(5-..>2)-T + H 


(45) 


Now we use Eq.(3) to get 
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lim o(^) = fi'o(^) = 21n 


(46) 



3 , 

H— In 
2 



27r2 

3 


+ 3.43 


This simple formula shows explicitly that there is no mass singularity in g{E), in complete 
agreement with the KLN theorem. 

An expression analogous to the RHS of Eq.(46) appeared in a footnote without proof 
in Ref. [9] It was given there as an asymptotic formula for G{Em) in the limit of large 
/3-decay endpoint energy, Em- The function G{Em), in turn, is the Sirlin’s function from 
Ref. [8] averaged over the f3 spectrum. The relationship between g{E) and G{Em) can be 
constructed from the dehnition of the latter: 


J,^-{Em-EfE^G{Em,E)dE 
j^-{Em-EYEME 
J,^-{Em-EfE^ {gGE)-2\n{^) 
ff-(Em - EfEME 


= 31n 



47r^ 


+ 7.1 


0 . 55 ) dE 


(47) 


The last line above coincides with the formula in Ref. [9] up to a small constant, which is 
unimportant for our discussion here. 

In Eq. (47), the limit rrie —0 is taken everywhere. The second line follows from the 
first line because the total beta decay rate must be the same in the limit cx) (the 

Erst line) and 0 (the second line). The constant contribution containing InM^ 

represents the “inner” radiative corrections, and it must be subtracted from g{E) if only 
the “outer” part of the correction is considered. Since the limit of large Em must formally 
coincide with the limit of vanishing rrie it is reassuring that both expressions have the same 
dependence on Em- 

Fig. 4 shows that for E > 1 MeV the functions g{E) and go{E) differ by no more than 
0.1%. This suggests, for example, that go{E) can be used in place of g{E) for the charged 
current deuterium disintegration reaction at SNO, where only the events with Eobs > 5 MeV 
were considered [3,4]. Incidentally, we have for large E [E is in units of MeV): 



3q; 

27rE2 


0.0035 

E2 


(48) 


in excellent agreement with Eq.(22). 


VII. CAPTURE REACTIONS 

Our analysis thus far has emphasized the universal features of radiative corrections for 
low-energy charged current processes in which the charged lepton appears in the final state. 


^The expression in Ref. [9] differes from Eq. (46) by a factor of two due to a different convention. 
We take this difference into account when comparing the two formulae 
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It is natural to ask whether the corrections relevant to capture reactions display similar 
universality properties. Examples of such reactions include: 


e + p —n + z/e 

P + P + e~ d + Pe 

^Be + e“ —> Li + z/g + ( 7 ) (49) 


The last two reactions produce pep and ^Be solar neutrinos, respectively. In the last reaction, 
the photon emission is due to ^Li nucleus decay from an excited state, to which ^Be decays 
about 10% of the time. Here, we show that an expression analogous to Eq. (15) applies 
for such capture reactions, with Fobs —^ Be and gbrem{Eobs) —^ 9b‘^^^iBe,Q), where Eg is 
the initial state charged lepton energy and Q is the Q-value for the reaction. The presence 
of g^°'^^{Ee,Q) implies that the correction factor in this case is not universal, though the 
non-universal dependence on Q can be computed in a straightforward manner. 

As far as radiative corrections are concerned, all reactions of the type shown in Eq.(49) 
can be treated in a unihed manner. Although at leading order the neutrinos are monoener- 
getic (for a hxed electron energy), bremsstrahlung from the initial state electron smears the 
spectrum to certain extent. Here, we will not study the spectral properties of the emitted 
neutrino but rather focus on the correction to the total emission rate. 

It is straightforward to see that the part of the radiative corrections due to exchange of a 
virtual photon is the same as the one contributing to g{Eobs)- Although the bremsstrahlung 
contribution is different, it can be easily obtained from Eq.(12) with help of crossing sym¬ 
metry and time-reversal invariance. Specihcally, one needs to change the sign of the photon 
4-momentum, which is equivalent to changing the sign of E^ and x in Eq.(12). To get the 
correction to the rate, we write down the analogue of Eq.(13) using the correct phase space. 
The result is: 

= rg 5 .(s., Q) x f -^)dxKdK , ( 50 ) 


where Q) is the tree-level capture rate for a given electron energy and Q-value of the 

transition. To study the correction to the total rate we can assume that the bremsstrahlung 
is never seen. This leads to the following expression: 




nl rEe+Q 

dx 


1-1 


2 In 


KdK ^^^t^ E^, -x) 


Ee + Q 

A 


mEe 


{Ee + qy 


1 - PiEe) 


+ C{P{Eq) 


{Ee + QP 1 f l+P{Ee) \ 
2AEI P{Ee) ""{l-PiEe)) 


llEe + 2Q 

^e 


mEe) 


( l+P{Ee) \ 

V-P{Ee)l 



(51) 


where the notation is the same as in the discussion of g{Eobs)- It is apparent that the part 
that contains the infrared divergence is the same as before. To get the total correction to the 
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rate we simply add the pieces that correspond to the virtual photon exchange (see Eq.(7) 
and Eq.( 8 )): 

gcaAEe, Q) = 9viE,r-^ + + g^^^\E,, Q) . (52) 

It is easy to show that the same formula applies to the positron capture case. 

To illustrate the dependence of gcapt on the Q-value of the transition we plot in Fig. 5 
the value of this function for a range of electron energies for three cases corresponding to 
the transitions in Eq.(49): 

Qp = Mp- Mn = -1.3 MeV , 

Qpep = 2Mp — Md = 0.93 MeV, ignoring pp kinetic energy , 

Qrse = Ml Be — M7Bi = 0.35 MeV, transition to the ground state of ^Li . 

(53) 

Note the turnover on the solid curve ((5p=-1.3 MeV). Generally, the turnover always appears 
if Q + me < 0. In such case the electron must have kinetic energy of at least TJ"*" = 
—Q — me for the capture to occur. The one-photon bremsstrahlung correction formally 
diverges at the threshold Eg -|- Q —^ 0, as evident from Eq.(51). This indicates that one- 
photon approximation breaks down, and contributions with arbitrary number of both real 
and virtual soft photons must be resumed to obtain a meaningful answer. In practical 
applications this nicety will rarely be of any importance, however, since the perturbation 
series breaks down only if 

a, ( me \ ^ ^ 

“ F I ^ oi Ee + Q < meC - (54) 

TT \Ee + Q/ 

In practice, the spectrum of the initial state electrons is always many orders of magnitude 
wider than the last term in Eq.(54). Since the leading order capture rate is suppressed by 
a factor of {Ee + Q)"^ the contribution from the immediate vicinity of the threshold to the 
average rate is negligible. 

In summary, the radiative corrections for capture reactions are not described by a uni¬ 
versal function. The initial state bremsstrahlung 7 - which is not generally detected - 
cannot radiate away more of the electron energy than required to make the reaction occur, 
thereby introducing a Q-dependence into the radiative corrections. Such considerations do 
not apply when the charged lepton appears in the hnal state, since in this case the min¬ 
imum detectable, observed energy is determined by the experimental conhguration rather 
the the Q-value of the reaction. The precise value of energy transferred from the lepton to 
the hadronic system does affect the energy-dependence of the total cross section, but this 
dependence is the same for both the leading-order and 0(a) contributions. Thus, apart 
from nuclear structure-dependent terms that are likely negligible for present purposes, the 
relative correction to the tree-level cross section for the reactions in Eq. (2) is described by 
a universal function. These features, then, have allowed us to formulate a unihed treatment 
of radiative corrections for neutrino reactions that we hope will help facilitate the analysis 
of future neutrino property studies. 
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APPENDIX A: CLOSED FORM FOR SOME INTEGRALS 

In Eq. (5) gb{E) is given in terms of two non-singular integrals over the electron energy. 
These integrals can be evaluated in closed form, with the results presented below. Our 
evaluation of the integrals is consistent with the results originally obtained in Refs. [22,23]. 
In practical calculations, numerical integration might, in fact, be more efficient than the 
evaluation of the exact expressions due to the complexity of the analytic formulas. However, 
closed form for the integrals may prove valuable if one is interested in studying the radiative 
corrections in various limiting cases (such as rrie —^ 0, etc.). For the first integral, we obtain: 


The second integral is more complicated: 



[^{E-x) In ( 

Jme y 


1 + /3{x) 
1 - /3{x) 
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where F{E) is defined in Eq. (5), and j3{E) and L{x) are defined in Eq. (6). 


(A2) 
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FIGURES 




FIG. 1. One-nucleon (a) and two-nucleon (b) contributions to the radiative corrections to CC 
neutrino-deuterium disintegration. The oval represents the initial state deuteron. 



FIG. 2. Exchanges of a low energy virtual photon for the processes shown in Eq.(2) in the 
one-nucleon approximation. 


P 



EIG. 3. The dominant bremsstrahlung graph. Ti and T 2 represent the initial and the final 
nuclear wavefunctions, respectively, q is the 4-momentum transfer from the leptonic to the hadronic 
part of the graph. In general, Ti ^2 can contain any number of nucleons. 
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FIG. 4. The exact one-loop radiative correction {a/TT)g{Eo^is) in % (solid line) for reactions in 
Eq.(2) and the same correction in the limit me —> 0 (dashed line). 



(MeV) 

FIG. 5. The electron energy dependence of the one-loop radiative corrections to the electron 
capture reactions in Eq. (49): p + e~ —> n -|- z^e (solid line), p-l-p + e~ d + i>e (dashed line), and 
^Be + e~ —Li + (dotted line). 
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